Liquid elasticity length and universal dynamic crossovers 
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We discuss two main universal dynamic crossovers in a liquid that correspond to relaxation times 
of 1 ps and 10 -7 — 10 -6 s. We introduce the concept of liquid elasticity length d e ;. At room 
temperature, d e i is several A in water and increases to 0.01 mm in honey and 1 mm in tar. We show 
that on temperature decrease, d e i crosses the fundamental lengths of the system, medium-range 
order d m and system size L. We discuss how d e i — d m and d e i — L correspond to the two dynamic 
crossovers. 



A conceptually simple phenomenon, freezing of liquid 
into glass, has turned out to be one of the most difficult 
problems in condensed matter physics, the problem of 
glass transition [l]. Analyzing the current state of the 
field, Dyre recently suggested that glass transition itself 
is not a big mystery: it universally takes in any liquid 
when its relaxation time r exceeds the time of experiment 
at the glass transition temperature T g 0. He proposed 
that the challenges lie above T g Q. 

If we consider the changes in molecular dynamics in a 
liquid on lowering the temperature, we find two dynamic 
crossovers. The first crossover is at high temperature at 
r « 1 ps, at which dynamics changes from exponential 
Debye relaxation to SER, q(t) oc exp(— (t/r)' 3 ), where q 
is a relaxing quantity and < (3 < 1. This crossover 
is universal, i.e. is seen in many systems 0, 0, H, 0]- 
SER describes a very sluggish dynamics: in the wide data 
range, it decays as a logarithm of time. As the tem- 
perature is lowered, we find another universal crossover, 
which takes place at lower temperatures that correspond 
to t = 1(T 7 - 1(T 6 s This crossover also marks 

the qualitative change in system's dynamics, and was at- 
tributed to the transition from the "liquid-like" to the 
"solid-like" behaviour. Note that although relaxation 
time at the second crossover is much larger as compared 
to the first one, it is still about 9-10 orders magnitude 
smaller that relaxation time at the glass transition. By 
definition, T g corresponds to r on the order of 10 3 s. 

One hopes to find useful insights into the problem of 
glass transition if the origin of dynamic crossovers could 
be rationalized. In this paper, we discuss this problem 
in terms of elastic stresses that a liquid supports in re- 
sponse to an external perturbation. We introduce the 
temperature-dependent liquid elasticity length, which is 
the range of elastic interaction in a liquid. We propose 
that that the first and second crossovers take place when 
this length becomes equal to the values of the medium- 
range length and system size, respectively. 

We first note that in a perfect crystal, there are two 
fundamental lengths, lattice constant a and system size 
L. In a disordered system like liquid, there is an addi- 
tional length d m , which corresponds to the medium-range 
order, and is defined by local packing. d m is on the order 



of 10 A, the characteristic size of decay of correlations in 
a disordered media. 

We now introduce the liquid elasticity length d e i- Un- 
like a solid, a liquid does not support static stresses. 
However, it supports stresses at high frequencies in a 
solid- like manner These frequencies correspond to 
times smaller than system relaxation time r. If d e i is the 
distance over which interactions in a liquid are elastic, 
d e i = cr , where c is the speed of propagation of elastic 
interactions, and is approximately equal to the speed of 
sound, c is on the order of c = a/ro, where To is the os- 
cillation period, or inverse of Debye frequency (to = 0.1 
ps) and a is the interatomic separation of about 1 A. We 
therefore find 



d e i = —a 
To 



(1) 



We note that d e i crosses, on lowering the tempera- 
ture, all three fundamental lengths in a system. Since 
t = ToCxp(Vb/fcT) at high temperature, where Vo is the 
high-temperature activation barrier, d e i = aexp(Vo/fc7 1 ). 
At high temperature, d e i is on the order of interatomic 
distance a. On lowering the temperature, it increases to 
d m . Because d m is about 10 A, we find from Eq. (1) that 
del = d m gives r of about 1 ps. This is the first dynamic 
crossover discussed above. On lowering the temperature 
even further, Eq. (1) shows that del increases to L. In 
liquid relaxation experiments, L is typically 1-10 mm. 
According to Eq. (1), d e i — L gives r = 10~ 7 — 10~ 6 , the 
second dynamic crossover. This points to the intrigu- 
ing possibility to describe the most important changes 
in system's dynamics using only its three fundamental 
lengths. 

To discuss why the two dynamic crossovers correspond 
to dei = d m and d e i = L, we first recall the previous dis- 
cussion how a liquid relaxes stress in response to an ex- 
ternal perturbation. Some time ago, Orowan introduced 
"condordant" local rearrangement events [13]. A con- 
cordant local rearrangement is accompanied by a strain 
agreeing in direction with the applied external stress, and 
reduces the energy and local stress (see Figure 1). A dis- 
cordant rearrangement, on the other hand, increases the 
energy and local stress. This has led to a general re- 
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FIG. 1: Orowan's illustration of a concordant local rearrange- 
ment. Solid and dashed lines around the shaded atoms cor- 
respond to initial and final positions of a rearrangement, re- 
spectively. Arrows show the direction of external stress. 



suit that stress relaxation by earlier concordant events 
leads to the increase of stress on later relaxing regions in 
a system ■ Goldstein applied the same argument to 
a viscous liquid 111 ]: consider a system under external 
stress which is counterbalanced by stresses supported by 
local regions. When a local rearrangement to a potential 
minimum, biased by the external stress, occurs (a con- 
cordant event) , this local region supports less stress after 
the event than before; therefore, other local regions in 
the system should support more stress after that event 
than before 



11] 



Lets consider stress relaxation in a liquid under ex- 
ternal perturbation (e.g., pressure or temperature). To 
illustrate the argument of Orowan and Goldstein above, 
we consider a system under stress and discuss the dynam- 
ics of pressure-induced local relaxation events (LREs), 
i.e those events that come in addition to thermal rear- 
rangements. In a simple model of parallel springs under 
external stress, a relaxation by a concordant event can 
be modeled as an appearance of a sp ring with a reduced 
spring constant. It is easy to show [121] , that after the 
introduction of n weaker springs, which corresponds to 
n concordant relaxation events in a real system, the in- 
crease of stress on a current unrelaxed spring Ap is pro- 
portional to the current number of events n, Ap oc n. We 
now invoke the previous discussion about the activation 
barrier V for a LRE in a real liquid. It has been ar gued 
that V is essentially the elastic energy 0, [ljj H3. Il5|. 
In the Shoving model, for example, V is given by the 
work of the elastic force needed to shove aside the sur- 
rounding liquid in order for a relaxation event to take 
place [15(| . Therefore the increase of activation barrier 
is proportional to the increase of elastic force due to the 
increase of pressure, Ap, on a current local rearranging 
region, and we find V oc Ap cx n. This provides the feed- 
forward interaction mechanism between LREs, in that 
activation barriers increase for later events. This mecha- 
nism is general, and does not depend on the choice of the 



spring model which we used here for illustration purposes 
only. 

It is important to discuss the condition under which the 
feed-forward interaction mechanism becomes operative 
on lowering the temperature. Let t s be the time needed 
for elastic interaction to propagate between LREs, d the 
distance between LREs, c the speed of sound, and To the 
oscillation period, or inverse of Debye frequency (to = 0.1 
ps). Because c = a/To, t s = Tod/a. Relaxation time t, 
which is also the time between two consecutive relax- 
ation events Q, increases on lowering the temperature. 
It is easy to see that at high enough temperature, when 
t = To, t s > t is always true because d/a > 1. In this 
case, local events relax independently of each other, be- 
cause at high temperature, the time between the events 
is shorter than the time needed for elastic interaction 
to propagate between them. Because events are inde- 
pendent, we obtain the expected high-temperature re- 
sult that relaxation is Debye (exponential) in time [l2| 
and Arrhenius in temperature. On cooling the system 
down, a certain temperature always gives the opposite 
condition: t s < t. When the time between local relax- 
ation events exceeds the time of propagation of elastic 
interaction between the events, the discussed above feed- 
forward interaction mechanism becomes operative. In 
other words, when t s < t, local relaxation events do not 
relax independently, but "feel" the presence of each other. 
We have recently showed how this picture gives rise to 
SER [12]. Note that d e i sets the maximal distance be- 
tween those LREs which are involved in the feed-forward 
interaction mechanism: t s = t gives d = cit/tq = d e i. 

We are now ready to discuss how the first dynamic 
crossover appears in this picture. Lets denote Ti the time 
of this crossover. As follows from the above discussion, 
Ti is equal to the time needed for elastic interaction to 
propagate between neighbouring LRE. Let d be the dis- 
tance between two adjacent LRE. Its value is on the order 
of 10 interatomic distances, or medium range order d m . 
Because t\ = d m /c and c = a/ To, we find 



Tl 



(2) 



We observe that t\ is the relaxation time of the system 
which corresponds to d e i = d m (see Eq. (1)). 

Because to = 0.1 ps and d m /a are roughly system- and 
temperature-independent, Eq. (2) predicts that t\ is a 
universal parameter, independent on temperature or sys- 
tem type. This is consistent with experimental findings 

% II H S 13- If d / a is on the order of 1° and T o = 0.1 
ps, we find from Eq. (2) that ti is about 1 ps, in good 
agreement with the experimental value of t\ in 1-2 ps 
range iiSiQ. 

To discuss the second dynamic crossover, we elabo- 
rate on how, on lowering the temperature, the increase 
of the radius of the feed-forward interaction increases ac- 
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FIG. 2: Illustration of the feed-forward interaction mecha- 
nism between local relaxation events. This mechanism oper- 
ates within the sphere of radius r — d e i/2, which increases 
according to Eqs. (1-2). Shaded and open circles represent 
local relaxing regions inside and outside, respectively, of the 
interaction sphere. 



tivation barriers for LRE. As discussed above, the max- 
imal distance between those LREs which are involved in 
the feed-forward interaction mechanism is equal to the 
elasticity length d el . According to Eq. (1), d el is on 
the order of interatomic distances at high temperature, 
but increases on lowering the temperature. The quali- 
tative argument, which is sufficient for this discussion, 
proceeds as follows, and is illustrated in Figure 2. Let N 
be the number of LREs, induced by external perturba- 
tion (i.e. those in addition to thermally-induced events) 
in the sphere of diameter d e i. As discussed above, remote 
concordant LREs result in the increase of stress, Ap, on 
a current local region. Because this increase is additive, 
Ap increases with N. Since N increases with d e i 7 Ap is 
a monotonously increasing function of c? e /. As discussed 
above, the activation barrier V oc Ap. Hence we find 
that V is a monotonously increasing function of d e \ . 

The second dynamic crossover originates in this pic- 
ture when, on lowering the temperature, d e i reaches the 
system size L. At high temperature, only small parts of 
the system are elastic, d e i <C L, as in the case of the first 
crossover when d e i = d m . When d e i = L on lowering 
the temperature, all LREs in the system are involved in 
the feed-forward interaction. Hence d e i = L marks the 
transition of the system from being partially to wholly 
elastic, and should manifest itself as a qualitative change 
in the liquid's dynamics. 

Another interesting prediction from this picture comes 
from the observation that at d > L, reducing the temper- 
ature further has a weak effect on V, because at this point 
the increase of V is related to the feed-forward interaction 
between temperature-induced events only, but not due to 
the increase V by way of the increase of d e i. As a result, 
the system is expected to show a crossover to more strong 
behaviour. Note that in this picture, the second crossover 
depends on system size. If, for example, the increase of 



V is described by the Vogel-Fulcher-Tamman law, the 
divergence at a finite temperature is avoided because as 
the system crosses over to more strong behaviour, the di- 
vergence is pushed to zero temperature. Only the infinite 
system does not have the second dynamic crossover, and 
has the divergence at a finite temperature. 

Experimentally, there is ample evidence for the dy- 
namic crossover in many systems. Most direct evidence 
comes from the sharp kink in the dielectric function 
Q. The crossover to the lower slope of relaxation time, 
with the effect that glass transition becomes retarded, 
is observed [l6j], in agreement with the prediction from 
our picture. Other experiments include NMR relaxation 
data, which detect a similar dynamic crossover [17|, the 
crossover in the relaxation of cage sizes in the positron an- 
nihilation experiments [l8| and changes of non-ergodicity 
parameter [19j. 

If t 2 is relaxation time at the second crossover, d e i = L 
gives (see Eq. (1)): 



L 

T 2 = —TQ 

a 



(3) 



Eq. (3) predicts that, similar to n, t 2 is a uni- 
versal parameter, independent on temperature or sys- 
tem type. This is consistent with experimental findings 
[1 SEE EI 111- 

Typical values of L used in the ex- 
periments are 1-10 mm, which is dictated mostly by the 
experimental setup. For example, smaller system sizes 
can be associated with surface effects, while larger sys- 
tem sizes can involve temperature gradients. Using this 
range of L and r = 0.1 ps, we find from Eq. (3) that 
t 2 = 10 -7 — 10~ 6 s, in good agreement with experimental 
values [1, Q. Note that at the second d e i = L crossover, 
1/t2 has the meaning of the eigenfrequency of the system. 

We note that relaxation time t 2 = 10~ 7 — 10~ 6 s cor- 
responds to viscosities rj of 10 3 — 10 4 Pa-s. These are 
much larger than room-temperature viscosities of famil- 
iar liquids like water, ethanol, olive oil or glycerol, for 
which rj = 10~ 3 — 10 1 Pa-s. For these liquids, Eq. (1) 
gives d e i in the range of 1-1000 Angstroms. Viscosity 
of honey (77 = 10 2 Pa-s) still falls behind the viscosities 
that correspond to relaxation time t 2 . Honey's elasticity 
length is d e i — 0.01 mm. Very viscous tar (77 = 10 4 Pa-s, 
d e i=l mm) has relaxation time close to t 2 and elastic- 
ity length comparable to the experimental length scale. 
Only extremely viscous pitch (77 = 10 s Pa-s, d e i=10 m) 
has relaxation time that exceeds t 2 and elasticity length 
in excess of the experimental length scale. Hence these 
examples illustrate that although the second crossover is 
way above the glass transition (see below) , it corresponds 
to quite high values of viscosity. They also give the feeling 
for the actual values of the introduced elasticity length 
in real liquids. 

In this picture, no dynamic crossover takes place at 
T„, consistent with the experimental observations. Ac- 
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cording to Eq. (1), r(T g ) — 10 3 s corresponds to d e i 
of 1000 kilometers, which explains why the two dynamic 
crossovers are seen way before T g is reached. In other 
words, the absence of a dynamic crossover in the vicinity 
of T g is due to the imbalance between our typical exper- 
imental times and sample sizes: at typical experimental 
time, the elasticity length is more than 8 orders of mag- 
nitude larger than the typical experimental length. In 
this context, we note that thermodynamic anomalies at 
T g are sensitive to cooling rates and observation times. 
Were the thermodynamic parameters to be measured at 
higher cooling rates and shorter experimental times, they 
would show the changes at temperatures that correspond 
to relaxation time r%. 

In summary, we have introduced liquid elasticity 
length, and showed that when it becomes equal to the 
medium-range length and system size, two universal dy- 
namic crossovers take place. The values of crossovers 
are fixed by a, d m and L, three fundamental lengths in 
a disordered system. From the physical point of view, 
the first crossover can be said to be more fundamental 
since it does not depend on system size L. Our finding 
represents a new interesting link between geometry and 
dynamics of the system. 
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